Sampling Distributions of Random Electromagnetic Fields in 
Mesoscopic or Dynamical Systems 
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We derive the sampling probability density function (pdf) of an ideal localized random electro- 
magnetic field, its amplitude and intensity in an electromagnetic environment that is quasi-statically 
time-varying statistically homogeneous or static statistically inhomogeneous. The results allow for 
the estimation of field statistics and confidence intervals when a single spatial or temporal stochastic 
process produces randomization of the field. Sampling distributions are particularly significant when 
the number of degrees of freedom v is relatively small (typically, v < 40), e.g., in mesoscopic systems 
when the sample set size TV is relatively small by choice or by force. Results for both coherent and 
incoherent detection techniques are derived, for Cartesian, planar and full- vectorial fields. We show 
that the functional form of the sampling pdf depends on whether the random variable is dimensioned 
(e.g., the sampled electric field proper) or is expressed in dimensionless standardized or normalized 
form (e.g., the sampled electric field divided by its sampled standard deviation). For dimensioned 
quantities, the electric field, its amplitude and intensity exhibit different types of Bessel K sam- 

(2) 

pling pdfs, which differ significantly from the asymptotic Gauss normal and X2 P ensemble pdfs 
when v is relatively small. By contrast, for the corresponding standardized quantities, Student t, 
Fisher-Snedecor F and root-_F sampling pdfs are obtained that exhibit heavier tails than comparable 
Bessel K pdfs. Statistical uncertainties obtained from classical small-sample theory for dimension- 
less quantities are shown to be overestimated compared to dimensioned quantities. Differences in 
the sampling pdfs arising from de-normalization versus de-standardization are obtained. 

PACS numbers: PACS: 02.50.-r, 03.50.De, 05.10.Gg, 06.30.Ka, 41.20.-q, 41.20.Jb, 42.25.Bs 



INTRODUCTION 



In the application of statistical methods to disordered 
and chaotic systems, particularly nonergodic and/or 
mesoscopic systems the role played by sampling 

distributions 0] is of fundamental importance because, 
in practice, data sets are necessarily of limited (and of- 
ten small) size N. Strictly, the Central Limit Theorem 
(CLT) is inapplicable to finite sample sets, a fortiori to 
small sets. Therefore, sampling distributions, rather than 
their underlying (parent) ensemble distributions, should 
always be used in any proper comparison between theo- 
retical and empirical probability distributions based on 
numerical data from practical measurement (experiment) 
or simulation (computation) , particularly when detecting 
or validating unexpected phenomena. The finiteness of 
the number of degrees of freedom for sample sets, can 
have a profound effect not only on the statistical mo- 
ments, but also on the functional form and shape of the 
probability density function (pdf) - particularly near its 
extremities - and on second- or higher-order statistics, 
such as auto- or cross-correlation functions and associ- 
ated spectral densities, etc. As a matter of fact, sam- 
pling pdfs draw on properties of both first- and second- 
order ensemble statistics: while the sampling distribu- 
tions are representations for local instantaneous sampled 



fields, the value of v as a sampling distribution parameter 
is governed by the correlation distance between points. 
This dependence on nonlocal properties of the field ap- 
plies to any physically limited realizable sampled region. 
For sample sets in which all N values are statistically 
independent (as we shall further assume throughout), 
v = N -I. 

Ensemble pdfs for ideal random 3-D electromagnetic 
(EM) fields are well known and have been amply inves- 
tigated in various physical applications, whether in un- 
bounded space Q or in the presence of an impedance 
boundary [9(. By contrast, their associated sampling pdfs 
have received little attention yet are pertinent and require 
characterization. Idealized random (stochastic) classical 
EM fields offer a paradigm for characterizing wave prop- 
agation or transport governed by dynamic multiple scat- 
tering (time- varying configuration or boundaries). Ex- 
amples include fields inside acoustic or EM mode-tuned 
or mode-stirred reverberation chambers (MT/MSRCs) 
- [Hj]; mesoscopic structures; random or turbulent me- 
dia; polarization and anisotropy in, e.g. , the cosmic mi- 
crowave background radiation [13|-[l5|; static or dy- 
namic optical diffusion using random phase screens; dif- 
fusing wave spectroscopy [I3|-[ll|; etc. 

For the estimation of statistics of the field, amplitude, 
intensity, energy density and power for such 'wave chaos', 
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the ensemble pdf of the field is usually adequate when 
the number of statistically independent partial contribu- 
tions and, hence, v approaches infinity. For example, for 
MT/MSRCs, the prerequisites are that the wavelength 
is very short relative to the physical dimensions of the 
cavity and that the observation time is long compared 
to the correlation time of the stirring process because, 
physically, a random EM field is governed by the spatial 
or temporal extent of the space- or time-limited process. 
On occasion, however, v is too small for these ensemble 
pdfs to be sufficiently accurate. In this case, the CLT or 
maximum-entropy principle may not be applicable. Nev- 
ertheless, in such low-dimensional cases, spatial and/or 
temporal averaging increases the value of v [23], which 
may have a substantial effect on the measured or per- 
ceived distribution of the field. 

A key issue addressed in this paper is that statistically 
homogeneous random fields are characterized by pdfs 
with distribution parameters (e.g., average, standard de- 
viation, number of degrees of freedom, etc.) whose values 
have to be estimated from the same set of data of the fluc- 
tuating field. As a result, these parameters themselves 
show sampling fluctuations because of the finite value of 
v. These variations give rise to bi- or multivariate fluctu- 
ations and increased uncertainty (i.e., wider confidence 
intervals) for the sampled field, compared to the corre- 
sponding ensemble distribution where these parameters 
are known exactly. 

The physical origin of the limitation of the value of v 
can be twofold. First, the potential (i.e., maximum at- 
tainable) number of statistically independent realizations, 
A max , may be practically restricted, even if an unlim- 
ited number of different states of the statistical ensemble 
were physically realized. For example, in a MT/MSRC, 
this case corresponds to the so-called undermoded regime 
(i.e., multi-mode operation but with less than, say, ten 
cavity modes being excited simultaneously) ; typically oc- 
curring at wavelengths where modal overlap is relatively 
small. As another example, in a relatively small sample 
of a random medium, the limitation on the value of v 
refers to the case of a relatively small loading fraction 
of inclusions. In this case, even the ensemble distribu- 
tion does not possess Gauss normal statistics. Secondly, 
A max may be unlimited but, for economical or other rea- 
sons, the sample size may have to be severely restricted 
(N <C A max ). In this case, while f ma x is potentially 
large, the value of v actually realized in a sample set is 
relatively small. 

In this paper, we investigate sampling pdfs of the com- 
plex value, amplitude and intensity of 1-D, 2-D and 3- 
D random EM fields that have ideal Gaussian ensem- 
ble probability distributions, for both the actual (non- 
standardized) observed EM quantities and for their stan- 
dardized forms. Here, the notions of standarized and nor- 
malized random variables refer to a random field quantity 
divided by its own (sample or ensemble) standard devi- 
ation or mean value, respectively. Coherent as well as 
incoherent detection techniques are considered in each 



case. However, specific nonstationary effects associated 
with the dynamics of multiple scattering, which are rel- 
evant to certain aspects of diffusing wave spectroscopy 
[l7l [Tl |. are not addressed here: each realized state of 
the system is considered in its quasi-stationary approxi- 
mation. The scenario investigated here is also relevant, 
in particular, to several practical applications in wire- 
less communications (e.g., mobile-to- fixed and mobile- 
to-mobile transmission), signal processing, wave propa- 
gation in turbulent media, modal noise in optical fibres 
under restricted-mode launch conditions [20|, etc. 



II. ELECTRIC OR MAGNETIC FIELD 



A. Coherent detection 



1. N on- standardized field 

Consider a modulated local analytic electric field E = 
E' —}E" received by an antenna (sensor) or scatterer im- 
mersed in a time- varying multi-scattering environment. 
A harmonic time dependence exp(jwi) is assumed and 
suppressed. If this field is made up of an arbitrarily large 
(theoretically infinite) number of fluctuating partial fields 
(i.e., modal or angular plane- wave spectral components) 
whose realization forms a random walk in the complex 
plane [23| , then, on account of the CLT (valid under very 
general but definite conditions [IS]), the associated con- 
ditional probability of E'v> coincides with the ensemble 
pdf, given by 



fE'O I (S E , 0) ,M E ,U) ) ( e I ( S £' w ' m E'<'1 ) ) 



2-7T S E ,(,) 



Here, M E ,o) and S E i(i) represent random variables in- 
duced by the sample mean and sample standard devi- 
ation of E'W , respectively, whereas lowercase symbols 
m E ,(,) and s E ,o) represent their corresponding sample 
values. If the mean and standard deviation of E'('> are 
known with certainty, with respective constant ensemble 
values li E /{i) and a E nn, i-e., 

f M E >(>) ( m s'<'> ) = *(%'f)-fe'(')) (2) 
fs E , (l) (s E >(>)) = 5(s E , w - a E ,o)) (3) 
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then the conditional pdf (p} coincides with the marginal 
pdf f E i(o{ e ) because, in general, 



SeK') f° r Gaussian E'^\S E ,(,) is a XpJV-i pdf, i.e., 
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with 



and because in (and only in) the case of a Gauss normal 
ensemble distribution of are the sample mean and 
sample standard deviation independent random variables 
0, i.e., 

fs E , {l) ,M B ,(,) ( S E>(>> > rn E ,(o ) = fs E , {l) (sb'C) )/m e , w (m E , ( o)- 
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From (dJ) and ©-©, we obtain 

(e'W 



The sampling pdf (J8j> is in its self-sufficient form [22j ]. 
i.e., it contains the standard deviation of S E i(o itself as 
(6) a distribution parameter. Alternatively, © can be re- 
expressed in terms of the ensemble standard deviation 
a E ,(,) of E'('\ because both statistics are related via 



exp 
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Here, Ojj/m = a E /V2 = y/p(r E tp, = v /C p/2a Ea where 
a = £, y or z for the three Cartesian components E a of 
E and where the value of p corresponds to the number of 
spatial dimensions in which E is being considered. 

Of the possible values of p, viz., 1, 2, or 3, its spe- 
cific value is governed by the polarization direction of 
the electric- field sensor and/or by any EM excitation or 
boundary conditions that may enforce a certain fixed (de- 
terministic) polarization state of the field. Thus, the 
Cartesian component of a vector field and the 3-D vec- 
tor field itself correspond to the cases p — 1 and p = 3, 
respectively. For a paraxially propagating polarized or 
unpolarized optical random field, p = 1 or p = 2, re- 
spectively. The case p = 2 refers to an unpolarized elec- 
tric field that is transverse to the local wavevector (i.e., 
randomly elliptically polarized field) and will be denoted 
by a subscript 't', whereas p — 1 corresponds to its po- 
larized detection or to a randomly modulated linearly 
polarized field, denoted by a subscript 'a'. These val- 
ues of p apply to local fields detected by an electrically 
small sensor, whose characteristic length is less than the 
spatial coherence length of the field, i.e., typically A/2 
in an unbounded medium. Larger values of p apply to 
electrically large detectors whose receiving cross-section 
(whether physical or as a synthetic aperture) is larger 
than A/2 in one or more dimensions. For example, for 
currents induced in a linear antenna of length L, we have 
p = max[l,Z/(A/2)]. 

If, instead of ([3]), Sjj'W exhibits random fluctuations 
or if its value is not known precisely, the sampling pdf of 



With the aid of [H, (3.471.9)], the sampling pdf of E'^ 
for /Zg;/(/) = in |5]) is obtained as a marginal pdf of the 
joint pdf f E , ( ,).s , M (e' (/) ,s B ,(,);7V), viz., 
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The pdf ifT^l) belongs to McKay's class of Bessel K dis- 
tributions [201 . These pdfs have also been obtained in 
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a variety of applied statisticalproblems, in electromag- 
netics and elsewhere; cf., e.g., [H[ and references in 
23] . It is further worth noting that the pdf (fT3|) has also 
been obtained in a different case, viz., as an ensemble pdf 
associated with an ad hoc Bessel K distributed field am- 
plitude (implying non-Gaussian/non-Rayleigh ensemble 
statistics) and derived via a Blanc-Lapierre transforma- 
tion Q • Results that are at least qualitatively compliant 
with those obtained from the pdf (fl"3|) have been observed 
recently in experiments involving correlated scattering 
@ . The pdf (ffU|) is also included in the class of Meijer G 
limit pdfs for complex fields in undermoded MT/MSRCs 
[23l ] . which in turn is a special case of more general Fox 
H distributions p9j . 

It is emphasized that the result (TIB"]) strictly holds for 
a physically ideal, viz., Gaussian random field if it is 
subjected to a finite-sized sampling (detection) process. 
Only in this case is a Gaussian marginal pdf f E ".')\s , (/) 
compatible with a non-delta-distributed S E i{i) . For the 
physical field itself, a Gaussian f E n>)\s , (() implies the 
number of degrees of freedom approaching infinity, hence 
resulting in a delta distribution for fs E , {l) j and vice versa. 

Fig. QJ, shows the pdf (ffli|) for selected values of N. 
Because of symmetry of f E /o) with respect to its mean 
value, only the positive half of the distribution is shown. 

If E'('> exhibits a deterministic, i.e., constant bias /i^'W 
as in ||7|), then 



JW>(e'W;iV) = 



0~E><> 



0~E«> 



K W -M S '(')| 
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(15) 



For nonnegligible "slow" fluctuations oi m E H') in CTJ, 
i.e., for a random bias (trend) of E'^'\ the results are 
easily generalized as follows. On account of linearity, the 
sampling variable M E ,(,) for Gaussian E'v> is also Gaus- 
sian with the same expected value but standard deviation 
cje'C) /V^V- Therefore, if the bias is itself random with 
Gauss normal distribution, i.e., 



exp 



/m„, w {rn E ,w ; N) 



A r ('» E /(,)-M E /(/)) 2 



2a A 



V 2 k/N a E , (l ) 
instead of ^j, then with ([S])-© we arrive at 
/ E , w (e'M;Af) = 

C E'C 
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For practical numerical integration of (jT7j) , the double- 
infinite range of m E ni) can be limited to, say, \jjl E i{o — 

5a E , (l )/VN,fJ-E'(') +5<j e «.')/VN]. 

For sufficiently large N, the field E'V* and its sample 
mean M E i{o are approximately independent (sharing N— 
1 degrees of freedom), whence E'^ — M E 'C) is then also 
Gaussian with zero mean and standard deviation o~ E , {l) ~ 

cr E /(/) yl + 1/N. Thus, ((T5|) with o E >(>) replaced by a* EI(l) 
serves as a first-order approximation to the pdf (fT7|) for 
large N. 



2. Standardized field 

For comparison of an experimental distribution against 
a theoretical sampling distribution, it is necessary to 
standardize the experimentally obtained E'^'K Then, 
rather than determining the pdf of E'v> , we require the 
pdf of the dimensionless variate 



X , W A E'^ - M E , W 



S E 'C) 



(18) 



i.e., for uncertain sample values of M E ,(i) and S E ,(,) that 
are unknown a priori. Their values are often to be esti- 
mated from the same limited sample data set of E'^ . For 
Gaussian E'('\ the M E ni) and S E H') ar e statistically in- 
dependent and their sampling distributions are Gaussian 
and XpN-i, respectively [30] - We further assume that 
the fluctuations of M E ni) are small compared to those of 
S E /(/) and, a fortiori, E'^'K 

It is shown in Sec. [XT] that for M B , W = ex- 

hibits a Student t distribution [3l| with pN — 1 degrees 
of freedom, i.e., 



f x ,o)(x' (,) ;N) = C xKl) 1 + 



A'Y 



-pN/2 



pN - 1 



(19) 



with X'W 



E'W/S EKI) for <£'W) 
='(') , 



0, where the sample 
as a ( 

ratio of the two sample values e'^'' and s E '(>) > and where 



value i'W = e'd/s E >(>) is estimated as a (dimensionless) 



C 



(^) V(pJV-i)t' 



r 



(20) 



The pdf (|19|) can be used to compare the empirical dis- 
tribution of the measured standardized complex field 
against its ideal theoretical sampling distribution, par- 
ticularly when N is relatively small (N ~ 40 or less). 
For multi-scattering contributions, the value of N can 
be associated with the number of scattering centers sur- 
rounding the receiver and/or the states during the mo- 
tion of the environment around a momentarily stationary 
receiver. 

Without loss of generality, we further consider a de- 
terministically unbiased field (fJ, E n_r) = 0). Fig. [T}d shows 
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the pdf of X'V) for selected values of N. Again, only 
the positive half of the distribution is shown. It can be 
seen that, compared to the asymptotic Gauss normal pdf, 
differences with f E no in Fig. [T£i occur mainly near the 
origin and in the tails. 



B. Incoherent detection 

Of significant practical importance, particular for mea- 
surements at optical wavelengths, is the case where the 
random field is sampled incoherently, i.e., by square-law 
detection of the electric or magnetic power, energy or in- 
tensity devoid of phase information. Yet one may wish to 
extract the sampling pdf of the complex-valued electric 
or magnetic field from such a scalar measurement. To 
this end, we borrow results derived in Sec. IIII Ai and re- 
strict further analysis to the case of a circular /e(e). (For 
an elliptic /^(e), data for a second linearly independent 
quadratic form of E' and E" are needed, e.g., E' 2 — E" 2 .) 
From the derivation in Sec. IIII Al it follows that the sam- 
pling pdf of J7'W = £/(') is (|2"o) after replacing p with 
p/2, because the ensemble pdf of U'^\S U ,(,) is \p whereas 
fiji(j) has a Xp_/v-i sampling pdf. With the variate trans- 
formation / B /(/)(e'W) = 2|e'W \f v , w (u'M = e'W 2 ), we ar- 
rive at 



| e /(/) | X ^[p(JV+l)-3] 



x%p(a-d-i] (vVM^-i)^J) . 

(21) 

For detection of a Cartesian component of intensity or 
power, i.e., p = 1, (12"T]) is equivalent with (Hi)) . 



III. FIELD INTENSITY (ENERGY DENSITY, 
POWER) 

For the field intensity U = U' + U" = \E\ 2 - as well as 
for the energy density or power, which are proportional 
to U - we can determine its pdf either on the basis of 
measured in-phase and/or quadrature components of the 
field (i.e., coherent detection, e.g., using a vector network 
analyzer) , or measured directly using a square-law detec- 
tor (i.e., incoherent detection, e.g., using a power meter, 
spectrum analyzer, field probe, etc.). We determine sam- 
pling pdfs for both cases. 



A. Incoherent detection 

1. Non-standardized intensity 

In many practical cases, the intensity or energy density 
is measured or perceived through a square-law detector or 



perception process. Compared to f E iw\s , w for coherent 
detection, the pertinent conditional pdf (cpdf) is now 
fu\S v - On account of the CLT, the sampling cpdf for 
U\S V is xl P , i-e., 



fu\Su (u\su) 



r(p) su \s v 



p-1 



exp ( -VP— ) ■ ( 22 ) 



The sampling pdf of Su, associated with an underlying 
circular Gauss normal E, is a xipN-i P^j i- e -' 



fs u (su;N) 




x exp —dpN — 



with 



pN 



(23) 



(24) 



Performing a similar calculation as in Sec. Ill A II and 
again using [25|, (3.471.9)] yields the sampling pdf of U 
as a marginal pdf of fu.Su ( u > s u I N). A Bessel K dis- 
tribution is again obtained, but compared to (| 1 3(1 it is 
now of a different type (i.e., the exponent of the power 
and the order of the Bessel function differ by a different 
amount), viz., 

/■+0O 

fu(u;N) = / fu\sM*u)fs (*»>N)**u (25) 
Jo 



Cu (u_ 



i p(JV+l)-| 



xK p{N _ x) y \ 2 x U/p[pN - - 




(26) 



with normalization constant Cjj given by 

2 Mp(n+i)-\] 



Cr 



T(p)T (pN 
x (pN - - 



2) 



i[p(iV+l)-i] 



(27) 



Again, (|26|) is one of McKay's Bessel K distributions. 
It is remarkable that for the case p = 1, the same type 
of Bessel K distribution (f26|) has been obtained in the 
context of sea echo for microwave radar but as an 
ensemble pdf of U a based on different starting assump- 
tions, viz., for a finite random walk in the complex plane 
that assumes a randomly fluctuating and large but finite 
number of steps (independent scattering contributions) 
that is distributed according to a negative binomial dis- 
tribution, as a discretization of a gamma distribution for 
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continuous stepping. For general values of p, the pdf 
([55]) was also obtained in [23| as a limit distribution for 
imperfect reverberation based on a Bayesian model for 
a physical process of omnidirectional scattering. Hence, 
the present derivation of ([26]) shows that Bessel K distri- 
butions are far more universal than previously thought, 
as they can arise under the much less restrictive condition 
of a mere small-sample effect for an underlying Gaussian 
field, as opposed to a need for any functional form of a 
priori distributions for a fluctuating is, whether chosen ad 
hoc or otherwise. However, any apparent departure from 
an underlying Gauss normal field distribution must be 
interpreted with due care and does not necessarily point 
to physical nonlinearity. The identification of Bessel K 
distributions as sampling distributions for small sample 
sets of received power is supported by measurements of 
the evolution of the distribution function of pulsed en- 
ergy in a reverberant cavity [32| . Such evolution is char- 
acterized by a steady growth in the number of multipath 
components and, hence, v. 

Figs. Hl, [3^,, and 0^ show the pdf of incoherently 
detected Cartesian (p = 1), planar (p — 2) and total 
(p = 3) field intensities, respectively, for selected values 
of N. The heavier tail and the sharper peak (mode) of 
the distribution at smaller values of N are characteristic 
features. Fig. shows the standard deviation of 
for these three cases. 



2. Standardized intensity 

For comparison with the empirical pdf, we consider the 
sampling pdf of the ratio of the two random variables U 
and Su, viz., 



W 



A U_ 

Su 



(28) 



The calculation of fw(w; N) is detailed in Sec. IA 21 The 
result is a Fisher-Snedecor F distribution 



, i , n C w ( w 
fw{w;N) = 

sw \sw 



p-1 



pN-h 



as a counterpart of (fl9|) . where 



Cw — 



nP/2 



T(pN 



■p) 



(pN 



T(pN 



(29) 



(30) 



The pdf ([29)1 for the incoherently detected and sample- 
standardized Cartesian (p — 1), planar (p = 2) and to- 
tal (p — 3) field intensities W a , Wt, and W is shown 
in Figs. H>, [3Jd, and 0)3, respectively, at selected values 
of TV. Comparison between the K distributed intensi- 
ties with their F distributed sample-standarized values 
in Figs. shows that, for a given value of N, the sam- 
pling pdf of W exhibits a more pronounced spread than 



for U, owing to the larger uncertainty of W caused by 
fluctuations of Sjj- Fig. [5Jd shows the standard deviation 
of ([29| for the three cases, indicating much larger stan- 
dard deviations for small N compared to those in Fig. 
[5^ for nonstandardized intensities. 



3. Non-standardized vs. non-normalized intensity 



In Sec. IIII A 11 fu(u; N) was obtained via the sample- 
standardized variate U / Sjj in the cpdf (122p . However, 
since fijj ^ 0, we could have equally used the sample- 
normalized variate U/Mu for the cpdf to reference the 
data in this case. To this end, instead of (|25| . we now 
use 



fu(u;N) = 

with m 

fu\Mu{u\mu) 
}Mu{mu;N) 



-Hoc 



fuwiu^mu) f Mu {m u ; N)dmu (31) 



Y(p) mu \mu J 



p-i 



(pN) 



( "mu 



T (pN) a Mu V o~M, 
mu 



exp 

pN-l 



-p (32) 

mu' 



x exp 



pN 



(33) 



and Hu = mjj, o~u — s/Nctmu: whence (13"TT) becomes 
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Ip(JV+l)-l 



xK 




p(N- 



with normalization constant C' v given by 

, A 2 p~M N + 1 ) a^p(n+i) 
Cjj = - 



T(p)T (pN) 



(34) 



(35) 



The form (|3"4"|) is readily re-expressed in terms of fiu by 
replacing o~u with Hu/y/p- Comparing ([34]) -([35 ]) to (f26|) - 
(|27[) shows that the non-normalized cpdf (|34|) is retrieved 
by replacing N—(2p)^ 1 in the non-standardized cpdf PU)) 
by N, i.e., resulting in a marginal increase of the num- 
ber of independent samples. In other words, the non- 
normalized cpdf is marginally closer to the asymptotic 
ensemble pdf (N — > +oo) than the non-standardized 
cpdf, resulting in slightly smaller uncertainties. This re- 
sult is made plausible by the fact that the uncertainty of 
the sampling mean value is smaller than for the sampling 
standard deviation. It is worth emphasizing that, even 
though fu(u; N) is for the (dimensioned) energy density 
U, the chosen route for arriving at this pdf - i.e., whether 
via intermediary standardization or normalization in the 
cpdf of U - has an effect on the number of degrees of 
freedom of the end result but not on the functional form. 
In summary, whenever pN 3> 1 (i.e., for most practical 
cases), standardization and normalization yield indistin- 
guishable final results. 
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4- Intensity of biased field 

Instead of (f22|) . the sampling cpdf of U\Su for the in- 
coherently detected intensity of a biased field is a gen- 
eralization of the so-called modified Nakagami-Rice m 
distribution, given in self-sufficient form by 



fu a \Sv a MsuJ 



x exp - ^l + 2k L 



- Ug + UgQ 



with 



|£aol 



Jo \ 2^/l + W L 



-y/U a0 U a 



(36) 



(37) 



where Cu a is a normalization constant. The general ex- 
pression makes allowance for the fact that E' a and 
may, in principle, have different pdfs (functionally 
and/or parametrically) , although in most cases the pdf 
of E a is circular, i.e., /#/ = /e"- The sample pdfs of 
Ut = Ux + Uy and J7 = U x + U y + U z follow similarly from 
two- and threefold convolutions of (jlfl)) , respectively. Fig. 
[5] shows the pdf for selected values of 2V. 

If U a is to be compared with measured data, then the 
sampling pdf can be similarly calculated from fx> and 
fx" as 



ff , n f° fx'SV^\N)f xl ,{^—^-N) A 

JwA w a) = W Q I ,- i ClX 

Jo V X V w a - X 

(41) 



A 



where W a = U a /S, 



u a - 



IV. FIELD AMPLITUDE 



where nu a is a sample value of vu a — Since 
o jj a — <Ju a -u a0 for any constant (deterministic) value 
u a0 , the pdf fsy (su a ; N) is still given by (|2"3"j) with p = 
1. Thus, 



Incoherent detection 



Non-standardized field amplitude 



fu a (u a ;N) 



C! 



vu a Jo 
x exp 



N- 



N 



x 



exp + 2k 



Io[2^T+2k L 



2) a Ur _ 
u a 



ds Ua . (38) 



Note that kjj a depends implicitly on su a : to first approx- 
imation, kjj a ~ SjJ a I 'fJ-Ua ■ 



B. Coherent detection 

In the case of estimating fu(u] N) from measurements 
of E, the sampling pdf takes a different form from that 
for incoherent detection, as we show next. Knowing 
the pdf of E a = E' a — jE'^, we can derive the pdf of 
U = \E\ 2 = p\E a \ 2 , where \E a \ 2 = E' a 2 + E? = U a is 
the intensity of a Cartesian field component, and simi- 
larly for the electric or magnetic power P oc U. Since 

/ ,(,) (u'c* ) oc / ,(,) (eL (,) = v/ uf)/ v/ we obtain 



The sampling pdf of A = y/ E' 2 + E" 2 follows in a 
manner similar to that for U. In this case, ^ISU has a 
X2p cpdf, given in self-sufficient form as 



fA\S A ( a \ s A) 



'_(np±il\ 
p V r (p) J 



r(p) s A 



2p-l 



x exp ■ 



V r(p) 



(42) 



which has the Rayleigh distribution as a special case for 
p = 1, while Sa has the X2pN-i distribution 



nN 1 ( r(pjV) . ^ 



2(pJV-l) 



(pN-i) T{pN-\)a SA 



x exp ■ 



P 2 \T(j>N-±) 




with as a an( i a A related via 



f Ua {u a ;N) = C L 



f u >(x;N)f u ,Au a -x;N)dx{39) 



f E ^;N)f E ,{^u7g-^;N) 
Cu a I 7= , ax 



(40) 



&S A = °A 



\ 



1 - 



1 / T(pN) 

P N-i \r(pN-$) 



(44) 



The sampling pdf of A is obtained as a marginal pdf of 
fA,s A ( a i S A] N) and is again a Bessel K distribution, but 
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of yet another type compared to (jT3|) and ([26)) , viz. 

r+oo 

f A (a;N)= / f A \ SA (a\s A )fs A (sA;N)dsA 
C A ( a 



<? A \v A 



xK p(N _ iy i I 2 



\ 




pN 



1\ a 



2y tr A 



(45) 

where is obtained, with the aid of [25|, (6.561.16)], as 



C A t 



r( P )r ( p n - 1) 




x pJV- 



i[p(JV+l)-l] 



|[p(JV+2)-l] 



(46) 



The pdf (|45|) for p = 1 has been obtained in [3| for the 
abovementioned scenario of a random walk with fluctuat- 
ing number of steps with negative binomial distribution. 
For general p, this pdf has been retrieved as a limit dis- 
tribution for imperfect reverberation in [23| . 

The pdf (|45|) for p = 1,2 and 3 is shown in Figs. [7^l, 
EK, and EH, respectively, for selected values of N. Fig. 
ITUk shows the corresponding standard deviations of (|2"9"|) . 
Compared to Fig. [5^, the increase of the standard de- 
viations with decreasing N is generally smaller and less 
dependent on dimensionality. 



2. Standarized field amplitude 

For comparison with the empirical pdf of a field am- 
plitude measured by an electric or magnetic field probe, 
we consider the sampling pdf of the ratio of the variates 
A and S A , viz., 



S A 

The calculation of /y(w; N) is detailed in Sec. 
the final result is shown to be 



(47) 
3l where 



f v (v;N) = C v v 



2p-l 



pN 



-\p(N+l) 



(48) 



with 




T{pN-\+p) 
T(pN - |)r(p) 



The sampling pdf (|48| can be referred to as a root-F 
distribution and constitutes a counterpart of Q19[) and 
(|29| . Fig. [TOb shows its standard deviation. 



3. Amplitude of biased field 

Unfortunately, unlike (|36|). the cpdf f Aa \s A ( a a\s Aa ) 
cannot be expressed in self-sufficient closed form, because 
a Aa for the Nakagami-Rice cpdf of A depends in a com- 
plicated manner on o E iu), Therefore, we do not fur- 
ther pursue this case. A simpler approach is to derive 
f A (a; N) via variate transformation of (|38|) . 



B. Coherent detection 

Using the variate transformation A = y/U, the pdf of 
the local field magnitude A = J A 2 X + A 2 y + A 2 Z = \E\ oc 



U follows. With A a = \JE' a 2 + E'f, (40|) yields 

fA a (a a ;N) 

= a a f Ua (u a = a 2 a ) 



= L Aa a a I — — ; = dx (50) 



X\ at, - x 



where C Aa is a normalization constant. The pdf for A a 
is shown for selected values of N in Fig. [TT] and for A in 
Fig. [TJ 

In fact, one may consider the amplitude of the sam- 
pling field itself as a marginal of the joint sampling 
pdf fE',E"(s',e";N), followed by variate transforma- 
tion to the amplitude A — + E" 2 and phase 
$ = tan -1 (E"/E'). Assuming that E 1 and E" are 
statistically independent (which, strictly, requires them 
to be normally distributed) so that fg/ E"( e 'i e"; N) = 
f E >(e';N)f E „(e";N), then 

f A (a;N) = ^V^ 1 f (| cos<^ - (a /a)\\ sin^" 1 



y/pN-1 



| a cos <p — ao 



KpN^ y/pN-1 



a sin < 



d(j) (51) 



(49) 



where ao = fiE' with fi E " = 0. 

V. EXTENSION TO NON-GAUSSIAN 
ENSEMBLE DISTRIBUTIONS OF THE FIELD: 
ITERATIVE BAYESIAN SCHEME 

From the theorem of total probability, it follows that 

/s_, w |E'f) « /b'M \s e , {0 Ss , w • (52) 
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For non-Gaussian E'^'\ for which the prior pdf fs E , (l) is 
more difficult to determine than in the exposition given 
before, ((52)) can be used in an iterative process by means 
of an update equation for fg /(/) , where the latter can be 

assigned an initial Xn-i distribution /l -* . This prior 
distribution, together with the non-Gaussian fEKO\s , w > 
then allows for calculating the posterior f s i(i) \e>o that 
serves as the "prior" pdf for the next iteration. Explicitly, 
denoting the ith iteration in this scheme by the super- 
script '(£)', the iteration process is specified by 

for i = 0,1 thus yielding f^™L,„ after n itera- 

j t> Je'W\s e , {i) 

tions, with similar relations for the intensity and ampli- 
tude. 

In general, (T52"|) - (f5"5]) can be used even with empirical 
pdfs /e'<.0\s , (;) ■ A similar equation and procedure can be 
used for determining the intensity or magnitude for inco- 
herent detection. This outline of a procedure is sketched 
here for sake of completeness and guidance. However, no 
explicit results are included because they require knowl- 
edge of fE'C)\s , (/ ) which must follow from a separate 
investigation, e.g., from experiment. Further study is 
needed with respect to the demonstration of convergence 
of ([55]) to a stable pdf. 

VI. CONFIDENCE INTERVALS 

95%-confidence interval boundaries from sampling pdfs 
of E, U, A as well as X, W, V and their components, are 
compared to corresponding boundaries for the ensemble 
pdfs in Figs. [T31 HH and [151 respectively. The Figures 
show that N ~ 10 or larger is needed in order for the 
sampling confidence interval boundaries to be close to 
the ensemble boundaries. For example, for an overmoded 
Fabry-Perot resonator, a resonator length of the order of 
five wavelengths or larger is needed to achieve this. 

VII. CONCLUSION 

In this paper, we studied sampling distributions for 
the analytic complex-valued field, the intensity (energy 



density, power) and the magnitude for Gaussian statisti- 
cally homogeneous random electromagnetic waves. The 
main results are 03|)-(III1), (fT7| . (fT9]) . (j26]) . (p9]) . (j45]) and 
(l4"5)) . A common feature is that lowering the number of 
degrees of freedom characterizing sampling distributions 
results in their tails becoming heavier and a consequent 
associated increase of widths of confidence intervals com- 
pared to the Gauss normal, \\ v an d Xiv ensemble distri- 
butions, respectively. Furthermore, it was shown that 
standardized quantities (i.e., the sampled random field, 
magnitude or intensity divided by its own sample stan- 
dard deviation, whereby the latter is itself considered as 
a random variable whose sample value is calculated from 
the sample data set itself) exhibit sampling pdfs that are 
characterized by a wider spread than for the case where 
that quantity has an a priori known (ensemble) standard 
deviation. The reason is that in the former case the ratio 
of random variables defines a bivariate sampling distribu- 
tion whose (univariate) marginal pdf is sought, whereas 
in the latter case the sampling pdf is univariate from the 
outset. The differences in sampling distributions arising 
from choosing normalization rather than standardization 
were shown to be negligible in all practical cases. 

While the focus in this paper was on the simplest of 
a priori chosen conditional pdfs for the local instanta- 
neous field in unbounded media, extension to some more 
general cases and boundary- value problems is straightfor- 
ward. For example, compound exponential distributions 
for anisotropic ideal random fields near a conducting or 
dielectric boundary |33| can be used instead, leading 
to the total (vectorial) intensity and amplitude of 3-D 
random fields near a planar isotropic surface by simple 
superposition. 
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APPENDIX A: SAMPLING DISTRIBUTIONS OF 
FIELD, INTENSITY AND AMPLITUDE 

1. Field 

Here we derive the one-dimensional sampling distribu- 
tion for the intensity \E\ 2 or energy density U ex \E\ 2 of 
a random statistically homogeneous Cartesian or vector 
electric (or, by extension, magnetic) field E. 

The pdf of an ideal random field E is circular with 
independent and identically distributed (i.i.d.) real (in- 

phase) and imaginary (quadrature) components E' = 

Re(-E) and E" = lm(E). As a result, these components 
can be studied in isolation from each other. The par- 
ent (ensemble) distribution of E is a central Gauss nor- 
mal distribution ((E) = 0). Therefore, the standardized 

TV-point sample variance Dn-i = (N—l)S 2 , (l) /cr 2 , ( , ) — 



T,f=i(K''- M E'0)) 2 /(T 2 EK , ) exhibits a Xn-1 P df > whilc 



y(') 



Di = [E'O — M S /(/)) 2 /<7^, (;) has a Xi pdf. Consequently, 
the ratio 



X 



, (/) a E'(') - M E , {I) 



D N . X /(N-1) 



(Al) 



has a Student t sampling distribution with N—l degrees 
of freedom, whence for E'^ itself, 



);iV 



r(f) 



N - 1 



='(') 



S E 'C) 



-N/2 

(A2) 



for N > 1, where e'^, m E ,(,) and s B /(/> are sample values. 
Compared to the sample value s E iw of the parent pdf, 
the sampling pdf exhibits an incre ased sampling standard 
deviation, viz., v / (N-l)/(N-3)s E , { , ) for N > 3. 

Typically, the Student t distribution arises in the char- 
acterization of the sample mean, i.e., for [(M E ,{,) — 

(E'^))/(a mi) /s/N)}/{S b k.) /a mi) ), instead of - 
^E'('))/ <t e'(')]/(Se'C) / a E'(o) as in (|Aip . Without lack of 
generality however, we further use the simplified defini- 
tion X'W = E'^/S E , ( >), i-e., M E , W = 0, because each 
sample set can always be centralized by its own sample 
mean value while maintaining its sampling pdf. 

The square of X'W, represented as Y'^/S Y , {I) = 

X'^ 2 = 2?'W /iS|,, w , exhibits a i 2 sampling distribution 
that follows from (IA2I1 as 



Wy' (/) ;iV) 



r(f) 



S Y Hi) 



,,'(') 



(JV-l)« y , w 



JV' 



-1/2 

(A3) 



which is a Fisher-Snedecor Fuv-i pdf of the ratio of 
the standard Xi variate D\ and the standard Xn-i var i~ 
ate Dat_i (both being statistically independent, because 
M E ,(,) and S^'to are independent for Gauss normal E'^), 
whereby each variate is divided by its corresponding 
number of degrees of freedom, i.e., 



t 2 = Fi, 



N-l — 



xi/i 



x^^/iN-iy 



(A4) 



2. Field intensity 

To find the sampling distribution of a squared Carte- 
sian or vectorial EM field with 2p i.i.d. components 
Yj (corresponding to a p-dimensional analytic complex- 



valued field), i.e., 
ZjDz, where 

- Y' + Y" " 



a 2 x , we consider the ratio 



Z : 



a Y >(>) 



E 



Y- Y/' 



ay cry y . 

2 — 1 V ' I— 1 



2p 



7, E 



(X, - m x f 



(A5) 



A' 



has a standard distribution, on account of the addi- 
tion theorem for 2p i.i.d. standard Xi variates, and 



A (2piV-l) S£, ( „ 



(A6) 



exhibits a standard x\ v n-\ distribution, with S 2 X , {I) = 

Y%=t (X' i (/) -M x , i , ) ) 2 /{2pN-l). Hence, their ratio G is 
a scaled Fisher-Snedecor F variate with (2p,2pN — 1) 
degrees of freedom: 



G 



2p 



Z/(2p) 



D z 2pN-lD z /(2pN-l) 
2p 



2pN - 1 ' 



-F, 



2p,2pN-l- 



(A7) 



With ct\e\ 2 = 2 v / pc r ^,(,), the ratio Z/Dz can be related 
to the standardized power or field intensity as 



S v 2JpS 2 E , 2^/p (2pN-l)S E ,/a E , 

2pN - 1 Z 



2VP D z' 

Thus, combining (|A7|) and |[A"8|) yields 

W = = y/p F 2p ,2pN-l- 



(A8) 



(A9) 



With (|A7[) . upon scaling the standard Fzp^pN-x distri- 
bution, the sampling pdf of G is 



fc{g\N) 



T{pN-\+p) gP- 1 



T(pN-±)T(p) (l + g f N -i 



(A10) 
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For the square of the Cartesian in-phase field component 
[Re(-E)] 2 or the quadrature component [Im(i?)] 2 (i.e., for 
P = 1/2), (|A10P reduces to (fA^l as expected. Finally, 
from (|A9j) , the sampling distribution oiW = U/ Sjj fol- 
lows as 



(2) 

valid for X2p distributions, from variate transformation 
of (|A11[) and noting that 



f w (w;N) = 



P 



r(pN-±+p) 

T(pN-l)T(p) (pN 

niP- 1 



w 

sw 



v 

sy 



s w 



sw 



, , dw = 2v-¥rdv, 



(A17) 



,/-</2 



(All) 



we obtain the sampling pdf of V at v — a/ sa as 
fv(v,N) = Cw 2 ?- 1 



This represents the pdf of the standardized received sam- 
pled power W at u>, i.e., for the ratio of sampled values 
u and Sjj- 

In the limit N — ► +oo, the exponent p(N + 1) — | 
in (|A11[) reduces to pN — i; the prefactor T(pN — i + 
p)/T( P N - i) = (pN - i) • . . . • (pN - | +p) for p > 1 
becomes (pN — \) p ; and sjj approaches the ensemble 
statistic aui whence 



pN 



-(pJV+p-|) 

(A18) 



which can be referred to as a root-i* 1 distribution, where 

p 



C 



A 



V 



( P N- 



fw(w;N)-+ fu(u) = 



P 



J>/2 



r(p)<7j/ \<tu 



p-i 



exp 



-Vp- 



<JJJ 



(A12) 



which is the ensemble \\ v limit pdf, as expected. The 
result agrees with a well-known limit theorem from 
probability theory Q stating that F m ^ n {x) converges to 
Xmimx) when n/m— + + 00, i.e., F q (m,n) — ► x 2 (m)/m for 
the corresponding quantiles. In the same limit iV^+oo, 
the asymptotic mean value of W is 



and the standard deviation for pN > 5/2 is 



(A13) 



aw = 



l2(2pN-l) 2 (2pN+2p-3) 



-P 



1. 



2p(2pN-3) 2 (2pN-5) 
Thus, we find that the coefficient of variation, i.e., 



(A14) 



W) 



l2pN+2p-3 
p(2pN-5) 



1 



— (A15) 



has its value reduced to that for the xip ensemble dis- 
tribution when N— > + oo, whereas for small N its value 
is substantially larger, indicating larger relative uncer- 
tainty. 
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For the standardized field mag nitude V = A/S_ 
V(S^/S W ) W, with 
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FIG. 1: (color on-line) Sampling probability density function 
for real or imaginary Cartesian component of electric field 
(p=l) at selected values of iV: (a) Bessel K sampling pdfs 
of E'J; [eq. (JT3J]; (b) Student t sampling pdfs of X' {,) = 
E'^/S ' /(o [eq. JT5J]. In both plots, the right tail becomes 
thinner for increasing values of N. 
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FIG. 2: (color on-line) Sampling probability density func- 
tion for intensity or energy density of Cartesian component 
of field (p=l) at selected values of N based on incoherent 
detection: (a) Bessel K sampling pdfs of U a [eq. (|26p ]; (b) 
Fisher-Snedecor F sampling pdfs of W a = U a /Su a [eq. ((29)) ]. 
In both plots, the right tail becomes thinner for increasing 
values of TV. 
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FIG. 3: (color on-line) Probability density function of inten- 
sity or energy density of planar field (p=2) at selected values 
of TV based on incoherent detection: (a) Bessel K sampling 
pdfs of Ut [eq. (I26|l ]; (b) Fisher-Snedecor F sampling pdfs of 
Wt = Ut/Sut [ e 1- l|29|l ]. In both plots, the right tail becomes 
thinner for increasing values of TV. 



10 



e 

X 




FIG. 4: (color on-line) Probability density function of inten- 
sity or energy density of total (vector) field (p~S) at selected 
values of N based on incoherent detection: (a) Bessel K sam- 
pling pdfs of U [eq. ([26}]; (b) Fisher-Snedecor F sampling 
pdfs of W = U/Su [eq. J29]}]. In both plots, the right tail 
becomes thinner for increasing values of N. 




FIG. 5: (color on-line) Sampling standard deviations as a 
function of N (a) for Bessel K sampling distributions (1261) of 
U a , Ut, and U; (b) for Fisher-Snedecor F sampling distribu- 
tions (Ull) of U a /Su a , and U t /S Ut , U/Su. 
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FIG. 6: (color on-line) Probability density function of in- 
tensity or energy density of Cartesian component of field 
(p=l) derived from coherent detection at selected values of 
N: Bessel K sampling pdf of U a [eqn. rt40jl ]. 
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FIG. 7: (color on-line) Sampling probability density func- 
tion of amplitude of Cartesian component of field (p — 1) 
at selected values of N based on incoherent detection: (a) 
Bessel K sampling pdf (J45J) for A a ; (b) root-f 1 pdf |g8]! for 
V a — A a /SA a - In both plots, the right tail becomes thinner 
for increasing values of N. 
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FIG. 8: (color on-line) Sampling probability density function 
of amplitude of planar field (p=2) at selected values of TV 
based on incoherent detection: (a) Bessel K pdf (|45[) for A t ; 
(b) root-f 1 pdf (0 for V t = A t /S At - In both plots, the right 
tail becomes thinner for increasing values of TV. 
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FIG. 9: (color on-line) Sampling probability density function 
of amplitude of total (vector) field (p=3) at selected values of 
N based on incoherent detection: (a) Bessel K pdf (|45(l for 
A; (b) root-f 1 pdf iggj for V = A/S A - In both plots, the right 
tail becomes thinner for increasing values of N. 
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FIG. 11: (color on-line) Sampling probability density func- 
tion of magnitude of Cartesian component of field (p=l) at 
selected values of N, derived for coherent detection. The right 
tail becomes thinner for increasing values of TV. 
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FIG. 12: (color on-line) Sampling probability density function 
of magnitude of total (vector) field (p=3) at selected values 
of N, derived for coherent detection. In both plots, the right 
tail becomes thinner for increasing values of N. 
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FIG. 13: (color on-line) Upper boundaries of 95%-confidence 
intervals for real or imaginary parts of E a , E t , or E. The 
lower boundary F _1 (0.025) = -_F _1 (0.975) is symmetric 
with respect to e = 0. 
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FIG. 14: (color on-line) Lower and upper boundaries of 95%- 
confidence intervals (a) for a ID Cartesian component U a 
(p—1), (b) for a 2D planar field Ut (p=2), and (c) for the 
3D vectorial U (p=3), normalized by the respective sampling 
standard deviations. For each line type, the lower and upper 
curves represent 2.5% and 97.5% percentiles, respectively. 
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FIG. 15: (color on-line) Lower and upper boundaries of 95%- 
confidence intervals (a) for a ID Cartesian component A a 
(p—1), (b) for a 2D planar field A t (p=2), and (c) for the 
3D vectorial A {p=3) , normalized by the respective sampling 
standard deviations. For each line type, the lower and upper 
curves represent 2.5% and 97.5% percentiles, respectively. 



